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Abstract 

There exists a large body of research on the lens of mammalian eye over the past several decades. The 
objective of the current work is to provide a link between the most recent computational models and some 
of the pioneering work in the 1970s and 80s. We introduce a general non-electro-neutral model to study the 
microcirculation in lens of eyes. It describes the steady state relationships among ion fluxes, water flow and 
electric field inside cells, and in the narrow extracellular spaces between cells in the lens. Using asymptotic 
analysis, we derive a simplified model based on physiological data and compare our results with those in 
the literature. We show that our simplified model can be reduced further to the first generation models 
while our full model is consistent with the most recent computational models. In addition, our simplified 
model captures the main features of the full model. Our results serve as a useful link intermediate between 
the computational models and the first generation analytical models. Simplified models of this sort may be 
particularly helpful as the roles of similar osmotic pumps of microcirculation are examined in other tissues 
with narrow extracellular spaces, like cardiac and skeletal muscle, liver, kidney, epithelia in general, and 
the narrow extracellular spaces of the central nervous system, the "brain". Simplified models may reveal 
the general functional plan of these systems before full computational models become feasible and specific. 


1. Introduction 

Biological systems require continual inputs of mass and energy to stay alive. They are open systems that 
require flow of matter, and specific chemicals, across their boundaries. Unicellular organisms can provide 
that flow by diffusion to and across cell membranes. Diffusion is not adequate over distances larger than 
a few cell diameters, i.e., larger than say 2 x 10 ~ 6 meters, to pick a number. For that reason, multicellular 
organisms cannot provide those flows to their cells by diffusion itself. Multicellular organisms depend on 
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convection to bring materials close enough to cells so diffusion to and across cell membranes can provide 
what the cell needs to live. 

The circulatory system of blood vessels-arteries, veins, and capillaries-provides the convection in almost 
all tissues. But there is one clear exception, the lens of the (mammalian) eye. The lens does not have blood 
vessels, presumably because even capillaries would so seriously interfere with transparency. The lens 
is large, much larger than the length scale on which diffusion itself is efficient. The lens must provide 
nutrients through another kind of convection, a microcirculation of water that moves nutrients into the 
lens and rinses wastes out of it. The microcirculation is in fact driven by the lens itself, without an external 
'pump'. The lens is itself an osmotic pump. 

The lens is an asymmetrical electrical syncytium in which all cells are electrically coupled one to an¬ 
other, with a narrow extracellular space between the cells (see Fig.l). The extracellular space is filled with 
ionic solution in free diffusion with the plasma outside cells. It may also contain specialized more or less 
immobile proteins and specialized polysaccharides, as well as containing obstructions formed by the con- 
nexin proteins themselves. The intracellular space behaves very much as a large single cell would, with 
the bio-ions of classical electro-physiology (Na + , K + , Cl) free to move without much resistance from cell 
to cell, and many solutes of significant size (say with diameter less than 1.5 nm) able to move as well. The 
intracellular media contains proteins particularly the crystallins responsible for the high refractive index of 
the lens. So the lens is an example of a bidomain tissue that has been studied in some detail, first in skeletal 
muscle, then in cardiac muscle, and syncytia in general. Electrical models of bidomain tissues have been 
developed and a general approach combining morphology, theory, and experiments has been applied in 
reference [1], showing how the lens could be studied in this tradition. 

A general approach to bidomain tissues was implemented [2] involving detailed measurements of 
morphology (best done with statistical sampling by stereo-logical methods [3]), impedance spectroscopy 
[4, 5, 6, 7, 8, 9, 10] using intracellular probes (micro-electrodes) that force current to flow across membranes 
to the extracellular baths [11, 1, 12, 13, 14, 15], electric field theory to develop models appropriate to the 
structure [16, 17, 18, 19, 20] analyzing the spectroscopic data with the field theory [21, 22] and checking 
that parameters change appropriately (i.e., estimates of membrane capacitance are constant) as extracel¬ 
lular solutions are changed in composition and concentration [20, 23]. This work was extended to deal 
with transport by Mathias and co-workers [24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40] and 
computational models of the water flow in the lens were later developed in some detail [25, 41, 42] and 
exploited with great success, reviewed in [43,44], also see [45, 46,47,48, 49]. 
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The original work on electrical models is cited here because it provides coherent support, involving a 
range of techniques and approaches, to the general view of syncytial tissues, used here and in later work. 
It also shows the range of approaches needed to establish a (then) new view of a tissue. 

Mathias [50, 51] realized that an asymmetrical electrical syncytium would produce convection, in par¬ 
ticular in the lens [31]: he and co-workers systematically investigated the flow of water, solutes, and current 
in the lens, which is (in our opinion) a model of interdisciplinary research, combining theory, simulation, 
and measurements of many types [24, 25, 26, 27, 28, 29, 30, 32, 33, 34, 35, 36, 37, 38, 39, 40]. Computational 
models of the water flow in the lens were later developed in great detail [25, 41, 42] and compared to the 
more analytical models. These models have been extensively tested and we are fortunate that comprehen¬ 
sive reviews have been written of great value to newcomers to the field, particularly [43, 44] as well as 
[45, 46,47,48,49,43, 41,44], 

Since the pioneering work on the models of lens microcirculation system proposed by Mathias et al. 
[51,21], numerous investigations have been carried out [52, 53, 23, 20, 19]. The microcirculation model has 
firstly relied on a combination of electrical resistance and current measurements and theoretical modeling 
[18, 54, 19]. More recently, in order to provide a better understanding the electric current flow and potential 
field, the detail structure of lens has been included [55, 41, 44, 43], describing the asymmetric biological 
properties of the lens and measurements of pressure have been made [47,28]. Different types of fluid flow 
[56,57] and transport properties of the ions have been introduced. Meanwhile, the lens model [55] has been 
extended to simulate age-related changes in lens physiology [58] and a variety of physiological processes 
[26, 59,60, 61]. Reviews of current studies on micro-circulation in lens are most helpful [30, 26, 33]. Despite 
this large body of experimental and theoretical work, it is not completely clear how they are related to each 
other. In particular, it is not clear how the latest computational models are related to the pioneering work, 
and how theoretical analysis is related to experimental findings. In this paper, we will provide such a link. 

Based on the microscale model for semipermeable membrane [62] and bidomain method [51], we con¬ 
struct a mathematical model to ensure that all interactions are included and treated consistently. Using 
asymptotic analysis, we derive a reduced model, which can be used to obtain most physiologically sig¬ 
nificant quantities except for the intracellular pressure. This simplified model can be further reduced to 
the model proposed by Mathias [51] with additional assumptions that Nernst potentials (that describe gra¬ 
dients of chemical potential of each ionic species) and conductance are are constant in space. However, 
we will show that neither the Nernst potentials nor the conductance are constant. On the contrary, they 
vary significantly from the interior to the surface of the lens Therefore, both of these quantities need to be 
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coupled as part of the solution. 

Our model also shows explicitly that the intracellular pressure is decoupled from the rest of the vari¬ 
ables. Evolution has chosen parameters so the intracellular pressure does not affect the other parameters 
of the lens in a significant way. They are robust to variations of intracellular pressure. The evolutionary 
advantage of this adaptation is not clear to us, but may be more obvious to other workers with a greater 
knowledge of clinical realities that show how the lens becomes diseased [63, 46, 48, 49]. Our simplified 
model suggests that all the quantities can be computed without knowing the intracellular pressure. On the 
other hand, we need to solve the full model to find the value of the intracellular pressure. Our model is also 
calibrated by experimental data and predicts the effects of gap junctions [28,47] described by a 'membrane' 
permeability jq n . 

Our new results extend but do not fundamentally change previous work on the lens. We strengthen the 
view that the lens provides an osmotic pump to maintain the microcirculation necessary to sustain a living 
lens, for the life of the animal. We imagine that similar osmotic pumps create microcirculation in other cells 
and tissues of the body. 

This paper is organized as following. The full model for micro-circulation of water and ions are pro¬ 
posed based on conservation laws in Section 2. In Section 3, we obtain the leading order model by identify¬ 
ing small parameters in the full model. Based on the boundary conditions and partial differential equation 
(PDE) analysis, a simplified version of the leading order model is proposed and compared with the existing 
models. The model calibration and simulation results are shown in Section 4. The conclusions and future 
work are given in Section 5. 

2. Mathematical model 

In this section, we present a 1-D spherical symmetric non-electro-neutral model for microcirculation of 
the lens with radius R by using the bidomain method [17]. The model deals with two types of flow: the 
circulation of water (hydrodynamics) and the circulation of ions (electrodynamics), generalizing previous 
bidomain models that deal only with electrodynamics. The model is mainly derived from laws of conser¬ 
vation of ions and water in the presence of membrane flow between intra and extracellular domains. We 
note that a similar approach may be useful in other tissues with narrow extracellular space, like the heart, 
cardiac muscle, and the central nervous system, including the cerebral cortex, 'the brain'. 

2.1. Water circulation 

We assume 
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Figure 1: Schematic diagram of lens. (A) The sphere of the lens with three landmarks: anterior pole (AP), posterior pole (PP) and 
equator (EQ); (B) The control volume in the bidomain model; (C) The micro structure of the lens: 1. intracellular region 2. extracellular 
region 3. cell membrane 4. gap junction (connexions); (D) Distribution of the gap junctions between the cell membrane at EQ or AP 
and PP; (E) A single gap junction which allows the water and ions flows. 
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• the loss of intracellular water is only through membranes flowing into the extracellular space, vice 
versa [17]; 

• the trans-membrane water flux is proportional to the intra / extra-cellular hydrostatic pressure and 
osmotic pressure differences, i.e. Starling's law, classically applied to capillaries, here applied to 
membranes [64]. In a system like non-ideal ionic solutions in which 'everything interacts with ev¬ 
erything else' [65, 66], this statement needs derivation as well as assertion. A complete, and rigorous 
derivation can be found in [62]; 

• in the rest of this paper, subscript / = in, ex denotes the intra-/extracellular space and superscript i — 
Na + , K + , Cl~ denotes the zth specie ion. 

Then we obtain the following system for intra and extracellular velocities in domain Q = [0, R] 

7 Jr {f MexUex ) 

— — ,L m ( Pex — Pin + 'Ymk]} r P(Qi n ~ Oex)) r (l a ) 

1 d 

— (r 2 (M ex ii ex + M in u in )) = 0, (lb) 

where zz; and P; are the velocity and pressure in the intracellular and extracellular space, respectively. And 
O; is the osmotic pressure with definition 

Oe X = E C L Om =E C m+ f'”, 

i i Vm 

where Cj is the concentration of zth specie ion in / space. y~ is the density of the permanent negative 
charged protein. In this paper, we assume the permanent negative charged protein is uniformly distributed 
within intracellular space with valence of z. Here At; is the ratio of intracellular area (/ = in) and extracel¬ 
lular area (1 — ex), A4 V is the membrane area per volume unit, j m is the intracellular membrane reflectance, 
L m is intracellular membrane hydraulic permeability, kg is Boltzmann constant and T is temperature. 

As we mentioned before, the intracellular space is a connected space, where water can flow from cell to 
cell through connexin proteins joining membranes of neighboring cells,, and the extracellular space is nar¬ 
row with a high tortuosity. The intracellular velocity depends on the gradients of hydrostatic pressure and 
osmotic pressure [62, 41, 51], and the extracellular velocity is determined by the gradients of hydrostatic 
pressure and electric potential [41, 67], 
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where ipj is the electric potential in the 1 space, T c is the tortuosity of extracellular region and }i is the 
viscosity of water, k e is introduced to describe the effect of electro osmotic flow, K; is the permeability of 
intracellular region (l — in) and extracellular region (/ = ex), respectively 

Thanks to Eq. 2, Eq.l can be treated as equation of hydraulic pressure. Due to the axis symmetry 
condition, homogeneous Neumann boundary conditions are used for pressure at r — 0. At the surface of 
lens r — R, we set the extracellular hydrostatic pressure to be zero and the intracellular velocity is consistent 
with Eq.2 


3 Ppr 3 Pi, 


dr dr 
Pex — 0 , 


-0, 


at r — 0, 
at r — R, 


-f 


( 3 ) 


= Ls ( Pin - 7s KbT (O in - Oex)) , at r = R, 


where q s is surface membrane reflectance and L s is surface membrane hydraulic permeability. 


2.2. Ion circulation 

With similar assumptions, the conservation of ion concentration yields the following ion flux system 


(4a) 
(4b) 

The ion flux in the intracellular region /? and ion flux in the extracellular region J l ex are defined as 


r 2 ^ r ( r2 -^ ex Jex) — ■M-vfmt 
^* r {r 2 (M ex J i ex + M in J\ n )) = 0, 
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(5a) 

(5b) 
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where D| is the diffusion coefficient of the zth specie ion in the / space. The Hodgkin-Huxley conductance 
formulation [68, 69] is used to describe the trans-membrane flux of ions across intracellular membrane and 
surface membrane 


jm ~ g Z i {tyin ‘Pex £ ) / 

Qi , s 

js — ~j ytpin ~ ‘Pex ~ Ej 


(6a) 

(6b) 


where E l = p^log ) 
zth specie ion. 


is the Nernst potential (an expression of the difference of chemical potential) of 


(A) 


(B) 



■ Na + channel 
K + channel 

■ Cl" channel 

■ Gap junctions 
■Na + /K + pump 


Figure 2: (A) Schematic diagram of ion circulation and the distributions of ion channels and pumps. The purple line represents the 
sodium circulation, the light green represents the potassium circulation and the brown line represents chloride circulation. The surface 
epithelial cells (dark black square) connect with the intracellular cells (light black hexagon) by the gap junctions (orange rectangle). 
The sodium and chlorine ion channels are located on the intracellular membranes, while the potassium ion channel and sodium- 
potassium ATP pumps are found only on the surface membrane. (B) Schematic diagram of water circulation. Trans-membrane water 
transport is through AQPO and AQP1 gap junctions. APQO gap junctions are located on the intracellular membranes, and AQP1 is on 
surface membrane. 


In Eq. 6, the intracellular ion conductance g l and surface ion conductance G l depends on the ion channel 
distribution on the membrane (see Fig. 2). Based on previous work [51, 17, 41], we assume that (1) only 
Na + and Cl - can leak between intracellular and extracellular through ion channels inside the lens and (2) 
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that there is no trans-membrane flux for K + between the extracellular and intracellular region, i.e. ]f n — 0. 

Similarly, homogeneous Neumann boundary conditions are used at r — 0. At r — R, the extracellular 
concentrations are fixed and Robin boundary conditions are used for intracellular concentrations due to 
the trans-membrane flux and pump. 


j lex = Jin = at r = 0 

l c 'ex = Cl, J\ n = f 5 + a', at r = R, 


where a 1 is active ion pump on the surface membrane. Here we only consider the sodium-potassium pump 
on the surface. The strength of the pump depends on the ion's concentration as in [29, 41], 
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Na 




( 8 ) 


where 
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C* 
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(9) 


Due to the capacitance of the cell membrane, assumptions of exact charge neutrality can easily lead 
to paradoxes because they oversimplify Maxwell's equations by leaving out altogether the essential role 
of charge. We use the analysis of [70] and thus introduce a linear correction term replacing the charge 
neutrality condition [51, 41], without introducing significant error (See also [71]), 
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(10a) 

(10b) 


where i] is the porosity of intracellular region and C m is capacitance per unit area. 

Multiplying each ion concentration equation in Eq. 4 with ez l respectively, summing up and using 
Eq.10, the sodium equations are replaced by the following equations 
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with boundary conditions 


dpex dp in 

dr dr 
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— G s p m — G'E 2 + fp, 

i 

where p in = {p in - p ex ) + \z\e& and p ex = {p ex - p in ) 


at r — 0, 
at r — R, 

at r — R, 


= Es i ' GS = E G ^ # = *E 
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In Eq. 11, we define the intracellular conductance <7j w and extracellular conductance cr ex as 


^ = (E D «( z! ) 2C «) /*Sn = (E D m( z ') 2G 


z 


It is obvious that system 11 might be derived using either Eq. 4 and Eq. 10. Therefore, we should drop 
either Eq. 4 or Eq. 10 when Eq. 11 is used. 

2.3. Non-dimensionalization 

Since lens circulation is driven by the sodium-potassium pump, it is natural to choose the characteristic 
velocity u* n by the pump strength n Na * 


— 


,Nn* 


o* 


( 12 ) 


where O* — 2 (C^ a + C ( f) is characteristic osmotic pressure. Using Eq.lb, we obtain the scale of u ex as 


* f—1 * 

u ex = d 0 U in , 


( 13 ) 
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where — H 1 - With the characteristic values for <p, P, C l , chosen as H C^ n + C5, we obtain the 

y vi j n i t K ex l c u u 

dimensionless system for lens problem as follows ( Detailed derivation is given in the appendix C.) 
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(141) 
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with homogeneous Neumann boundary conditions at r = 0 and following boundary conditions at r — 1 


where 


P ex — 0, 
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(15a) 

(15b) 

(15c) 

(15d) 

(15e) 

(15f) 


3. Simplified model 

The full model given by system 14 with boundary condition 15 is a coupled nonlinear system. In this 
section, we present a simplified version of the full model which captures the main features of the lens 
circulation. We first obtain the leading order model by identify the small parameters. And then by using 
boundary conditions and theoretical analysis, the leading order model with is further simplified as only 
one PDE with serial algebra equations. 

According to those dimensionless parameters presented in the appendix B, we identify the scale of the 
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parameters as follows 

{44}cO(e), {^ 0 /^ 3 } C 0(e 2 ), 

{*2/ *4/ *5/ ^6/ <4} C o(e 2 ). 

If we denote £9 = Dp* — Df and *10 — Dp* — D ; Nfl , 1 — in, ex, it yields 

<S 9 = 0(e 2 ), *10 = 0 ( 6 ). 


(16) 


(17) 


3.1. A priori estimation 

In this section, we provide the priori estimation of the /p* as follows. By using the homogeneous Neu¬ 
mann boundary condition at r — 0 and Eq. 141 yields 

dr*™ = oP ( Pe inPinUin + 6 9 - c f n + 5 w -C? n a ^ 
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From Eq. 18, since Pe ul — O(e) and order of (5s, * 9 , *10 in Eqs. 16-17, we obtain that 
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Meanwhile, from Eq. 14b we can have 

and in the Eq.l4e, we know 



With Eqs. 20-21 and p" is constants, we obtain 

1 v in 
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Furthermore, using Eq. 14d and boundary conditions for Cp r * in Eq. 15 yieldds 
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From the experimental setting of lens [51, 55, 41], we assume that 


C 0 Nfl + Cf - TT = °( £ )- (24) 

Therefore, 

C% = 0(e). (25) 

In all, we claim that 

I hi - Pe in C£ «m - D ; 0 ^ Cg + Z° Cg j r (pin ) ^ 

— 0(e 2 ). 

By dropping the terms involving these small parameters, the leading order of water circulation system 
14a-14d is as follows. 
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where the superscript '0' denotes the leading order approximation. From Eq. 27, we deduce 0® x = O? are 
constants, and the intracellular and extracellular flow are counterflow. And the total charge in the leading 
order systems are neutral 


E 2,c h+4r = o, 

l v in 
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(28b) 


i 


Combining constant osmotic pressure and charge neutrality yields 





Cl 0 Cl 0 

which means C- ' and C ex ' are constants and 
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And the leading order of potassium and chloride concentrations satisfy 
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where / ; ! '° = Pe/C|'°u° — DJ + zl C° with i — K, Cl and / — zn, ex. 

For the electric potential, using the homogeneous Neumann boundary condition at r = 0 and Eqs. 29a 
-30,141 yields 


drQin ~ } in ( P^inPin u in + h dr^fn + ^10 drtfn') 


(32) 
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At the same time, based on the intracellular equation of potassium Eq. 14j the homogeneous Neumann 
boundary condition at r — 0 and Eqs. 29a -30, we have 
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Substituting Eq. 32 into Eq. 33 yields 
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where we used that fact that p,, r — 


Pex — o(e 2 ), S 9 — 0(e 2 ) and = -‘fojfo + 0(e 2 ), 1 e {in, ex}. Since 


15 







Pei n — 0(e), and S$ — O(e), in Eq. 34, we claim 
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Combining Eqs. 32 and 35 yields the leading order approximation of intracellular potential 
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where cr° = £ D^z’) 2 ^, < = D D^z*') 2 ^. 

Similarly, the leading order approximation of extracellular potential is 
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To summarize, the leading order approximation of system 14-15 is given by, in domain O — [0,1] 
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with boundary conditions at r — 1 
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In the following, we will further simplify Eqs. 38-39 and obtain the relationships between tp® x and other 
leading order variables by using assumptions concerning the boundary conditions. 


3.2. Relation between <pf n and (p® x 

Combining Eqs. 38a, 38d and 38k, and integrating with respect to r yields the relation between (p® and 


Tex as 
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where we used the boundary conditions <p® x (1) — P,? Y (1) = 0. 

3.3. Relation between P® x and <p® x 

By the homogeneous Neumann boundary condition on r — 0 and Eq. 38j, we have 
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Cl 0 

By Eq. 29b, we can divide Eq. 41 by C ex ' on both sides, we get 
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= 0 . 


Base on the charge neutrality Eq. 28, constant osmotic pressure Eq. 29a and parameters in Appendix B, 
we denotes 
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Then combining the Eqs. 36 and Pe !n = 0(e), Eq. 42 yields the following equation by omitting the 
higher order terms 
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Finally, by using the boundary condition, we have the relation between extracellular pressure and electric 
potential as 

D ex - Pe exS 1 ^ 


P U = 

x ex 


Pe P 


ex- 


3.4. Expression of E Na 

Based on potassium equation and relation in Eqs. 40 and 45, we have expression for Cf n and Cf x as 
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where 


Based on Eq. 28, we can get 
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3.5. Extracellular electric potential system 
By Eqs. 40 and 45, we have (p m as 


\ a in a in J 

The value (pf ( 1) is determined by the boundary condition of (p® n in Eq. 39, where 
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To summarize, we obtained the simplified model of system 38-39 as follows 
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(51k) 
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Remark 3.1. Under the same assumptions in [51], for example, uniform diffusion constants for all ions, constant 
Nernst potential, our simplified model system 51 recovers the model proposed by Mathias. The main contribution here 
is that we remove the assumptions that Nernst potentials and effective conductance shoidd be constants. By using 
the relationships between ions concentrations and external potential, zve obtain the space dependent Nernst potential 
zvhich yields a much better approximation to the full model (see Fig. 4). 

4. Results and discussion 

In this section, we present numerical simulations using both the full and simplified models. Finite 
Volume Method [70] is used in order to preserve mass conservation of ions. The convex iteration [72] is 
employed to solve the nonlinear coupled system. The numerical algorithm is implemented in Matlab. 

4.1. Model calibration: membrane conductance effects intracellular hydrostatic pressure 

In this section, we first calibrate the full model by the comparing with the experimental data to study 
effect of connexin to intracellular hydrostatic pressure. 

Intracellular hydrostatic pressure is an important physiological quantity [63]. In the paper [28, 47], the 
authors showed the connexin (gap junction) conductance play an important role in the microcirculation of 
lens. It is said that if the intracellular conductance jf 3 - in lenses is approximately doubled, the hydrostatic 
pressure gradient in the lenses should become approximately half of the original one. In this section, we 
calibrate our model. We choose a value of the intracellular conductance (—) that correctly calculates the 

Pin J 

experimental results in the [28, 47]. 

In Figure 3 (A), the value k® = 4.6830 x 10 ~ 20 /nr (black line) yields a good approximation to exper¬ 
imental data (black makers). When the conductivity of the connexins is doubled, to parameter value K ln 
to be 2 k® (in the lens of mice Cx46 KI lens) as in the experiments [28,47], where doubled the conductivity 
of the connexins by using Cx46 KI mice lens, our model (black dot) can also match the experimental data 
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(red markers): the intracellular hydrostatic pressure drops to half. This result shows that our full model 
can correctly predict the effect of permeability of membrane on hydrostatic pressure. 

Interestingly, panels (B)-(D) it shows that other intracellular quantities and extracellular ones (appendix) 
are insensitive to increases in the permeability by a factor of twenty, even to 20 k ® . The reason for this can 
be explained by using our simplified the system 51. If the variation of intracellular conductance still keep 
the &2 to be a small quantity in the dimensionless system 14, our simplified model will be still valid. In the 
simplified model. All the quantities except intracellular hydrostatic pressure are related to the extracellular 
electric potential. However, the extracellular electric potential will not be effected by the change of the 
intracellular conductance, since Eq. 51a not involves intracellular conductance. 



Figure 3: Comparison between different jc;„. The experimental data of dog, rabbit, rat come from paper [47], Mice and Cx46 KI mice 
come from paper [28]. According to paper [28], the Cx46 KI mice lens has twice the number density of lens gap junction channels 
compared to mice. The parameter jtF = 4.6830 x 10~ 20 / m 2 and radius is written in dimensionless units for different species 


4.2. Full model vs simplified model 

In this section, we compare the full model 14-15 with the simplified model 51 and Mathias model in 
[51]. The numerical results of full model (Black lines) in Figure 4 (A-C) suggest that the variations in¬ 
tracellular electric potential, extracellular conductance and Nemst potential of Cl are rather small. The 
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assumption of constant values for those variables ( potential, extracellular conductance, and Nernst, i.e. 
chemical potential of Cl - ) in the Mathias's model (shown as red dash-dot lines) is reasonable. However, 
the Nernst potentials of sodium and potassium (Figure 4 (D-E)) have large variations, because of the effect 
of Sodium-potassium pump. Our simplified model (black dash lines) describes these variations with small 
errors. The comparisons for extracellular pressure, velocity, and potential (Figure 4 (F-H)) confirm that our 
simplified model yields good approximations to the full model. 


5. Conclusion 

In this paper, we propose a bidomain model to study the microcirculation of lens. We include a capac¬ 
itor in the representation of the membrane and so our model is consistent with classical electrodynamics. 
Consistency produces a linear correction term in the classical charge neutrality equation. This full model is 
calibrated by comparing with the experiment studying effect of connexin on hydrostatic pressure. It shows 
that only by changing intracellular membrane conductance (strength of connexion), our model could match 
the two experimental results with different connexin very well. Our model is capable of making prediction 
to the circulation of lens. Furthermore, the numerical simulations show that the velocity, potential, osmotic 
pressure in the intra and extra cellular are not sensitive to increasing conductance. 

Based on the asymptotic analysis, we proposed a simplified model, which allows us to obtain a deep 
understanding of the physical process without making unrealistic assumptions. Our results showed that 
the simplified model is a good approximation of the full model where Nernst potentials and conductivity 
vary significantly inside the lens. 

Our model allows calculation of variables that determine the role and life of the lens as an organ. Partic¬ 
ularly important are the factors that determine the transparency of the lens, since that is the main function 
of the organ. The dependence of the size of the extracellular space, and thus the pressure in the extracellular 
and intracellular spaces and the difference between those two, is likely to be an important determinant of 
transparency. One imagines that swelling of the extracellular space will scatter light, particularly because 
the swelling is likely to be irregular (in a way our model does not yet capture). Changes in the Osmolarity 
(i.e., activity of water estimated by the total concentration of solutes) is likely to be important as well. 

This hydrodynamic bidomain model can point the way to dealing with other cells, tissues, and organs 
in which current flow, water flow, and cell volume changes are important. These include the kidney, the 
central nervous system (where the narrow extracellular space poses many of the biological problems facing 
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Figure 4: Comparison of electro-neutral and simplified and Mathias's model in [51]. 
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the lens), the t-tubular system of skeletal and much cardiac muscle and so on. We show that a mathemat¬ 
ically well defined model can deal with the reality of biological structure and its complex distribution of 
channels, etc. 

Conservation laws applied to simplified structures are enough to provide quite useful results, as they 
were in three dimensional electrical problems of cells of various geometries [16] and syncytia [3]-[23] . The 
exact results are analyzed with perturbation methods, described in general in [73] and these methods allow 
dramatic simplifications without introducing large or even significant errors. It is as if evolution chose sys¬ 
tems in which parameters and structures allow simple results, in which parameters can control biological 
function robustly. 

Of course, we only point the way. Additional compartments and additional structural complexity will 
surely be needed to deal with the workings of evolution. But these can be handled in a mathematically 
defined way, yielding approximate results with clear physical and biological interpretation. Combining 
the multi-domain model and membrane potential dependent conductance, one can model depolarization 
induced by extra potassium in lens [53, 55] and cortical spreading depression (CSD) problem [74, 75, 71]. 
The ultimate goals will be (i) to provide as much precision in the mathematics and physics as we can, start¬ 
ing from first principles [62]; (ii) to provide a general basis for treatments of convection in other tissues that 
involve microcirculation. Computational models of these are not in hand, and may be hard to construct, 
since so little is know of those systems compared to the lens. With what we have learned here, we hope a 
general mathematical approach and model of the type we present here may be constructed and helpful in 
other systems with narrow extracellular spaces that are likely to need microcirculation to augment diffu¬ 
sion, like cardiac and skeletal muscle, kidney, liver, epithelia, and the extracellular space of the brain. 
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Appendix A. Model Parameters 


Parameters 

Mathias [51] 

Malcolm [55] 

Parameters 

Mathias [51] 

Malcolm [55] 

R 

1.6 x 10~ 3 m 

1.6 x 1(T 3 m 

Pin 

3.75 x 10- 13 m/(Pa s) 

1.34 x 10 - l3 m/(Pa-s) 

Ai„/V in 

78 mM 

78 mM 

L s 

3.75 x 10~ 13 m/ (Pa-s) 

8.89 x 10 - 13 m/(Pa s) 

C Na 

107 mM 

107 mM 

Mi„ 

0.988 

0.99 

C K 

''■'O 

3 mM 

3 mM 

Men 

0.012 

0.01 

c„, 

- 

1 x ltr 2 F/m 2 

M„ 

6 x 10 5 /m 

5 x 10 5 /m 

D Na 

ex 

- 

1.39 x 1CT 9 m 2 /s 

T 

- 

310 K 

d k 

‘-'ex 

- 

2.04 x 10~ 9 m 2 /s 

k e 

1.72 x UT 8 m 2 /(V ■ s) 

1.45 x 10~ 8 m 2 /(V-s) 

D a 

‘-'ex 

- 

2.12 x 10~ 9 m 2 /s 

kB 

1.38 x Ur 23 I/K 

1.38 x 10- 23 J/K 

D Na 

- 

1.39 x 10~ n m 2 /s 

K K i 

- 

1.6154 mM 

D i 

- 

2.04 x 10~ n m 2 /s 

Kk2 

- 

0.1657 mM 

Dg 

- 

2.12 x 10- 11 m 1 /s 

KtsJal,Na2 

- 

2.3393 mM 

e 

1.6 X 1(T 19 A ■ s 

1.6 x 10~ 19 A ■ s 

>7 

0.988 

0.99 

g Na 

2.2 x 1(T 3 S/m 2 

2.2 x 10~ 3 S/m 2 

Kex 

1.141 x 10~ 16 m 1 

1.33 x 10~ 16 m 2 

g CI 

2.2 X 1(T 3 S/m 2 

2.2 x 10~ 3 S/m 2 

Kin 

- 

9.366 x 10~ 19 m 2 

g k 

2.1 S/m 2 

2.1 S/m 2 

7m,s 

1 

1 

h 

2.3 x 1( r 2 A/m 2 

- 

T c 

0.16 

0.16 

Imaxl 

- 

0.478 A/m 2 

i 1 

7 x 10 -4 Pa ■ s 

7 x 10 -4 Pa ■ s 

Imax 2 

- 

0.065 A/m 2 

2 

-1.5 

-1.5 
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Appendix B. Dimensionless Parameters and Scales 


The following dimensionless parameters' value and scales calculation based on values in [55] 


Scales / Parameters 

Value 

Parameters 

Value 

a N «* 

6.9 X 10 - 2 A/m 7 

& 

H 

aa 

1 

99 

C* 

110 mM 

X k e T C k B T 

1 

1.2031 x 10- 1 

O* 

220 mM 

„ uRu* 

— * in 

2 K in y m k B TO* 

6.861 x 10- 3 

P * 

16.937 KPa 

x _ P* 

3 7mk B TO* 

2.9894 x 10~ 2 

U in 

3.2506 nm/s 

Si 

3.5323 x 10~ 5 



M in«t 




RM v Lm7mk B TO* 


Kx 

3.2181 }im/s 

^ L sls k B TO‘ 

4.3022 x 10~ 3 

F 

26.7 mV 

X MvCmk B T 

1.2745 x 10- 5 

D ex 

3.392 x 10 10 m 2 /s 

<- MvCmk B T 

52 - e 2 C*(l— i/} 

1.2617 x 10~ 3 

Vn 

2.12 x 10~ n m 2 /s 

c _ MexDg X 

8 M mDl 

1.6162 x lO^ 1 



Df'-Df 


Pe ex 

1.5180 

v 

3.77 X 10 2 


2.4533 x lO^ 1 



Pe„ 

* 0 = 39 - 

cP'° 

= TOT 

'-ex 

3.443 x 10 1 


d n« 

in,ex 

0.6557 

12.5 

110 

fb 

in,ex 

0.9623 

Po 

117 

110 

m 

in,ex 

1 

M<" 

3.3859 x 10- 1 

R s 

4.00 x 10- 1 

M “ 

2.095 


The 8 can be find in the following equations. 


4: in e^. [13], 

<5 7 : in e^.[14/]. 

8 \ : in eq.[14a]. 

: in eq.[14h], 

S 2 : ineq.[ 14b], 

89 : in eq. [18], 

83 : ineq.[ 146], 

^ 10 : in [18], 

8 ±\ in eq. [14c], 

8 n : in eq.[ 43], 

: in B.C. beiorc [14], 


: in e^.[14e], 
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Appendix C. Non-dimensionalization 

In this section, we derive the dimensionless model based on the lens, which has been widely stud¬ 
ied. The major ions we considering here are sodium ( Na + ), potassium ( K + ) and chloride (C/ _ ) and the 
sodium-potassium pump which distributed on the surface of the lens. Although we restrict ourselves in 
this particular problem, the following procedure can be applied in a wide range of practical problems in 
biological syncytia. 

Appendix C.l. Water circulation 

In the following, we assume the typical length of lens is R. The fluid system is driven by the osmotic 
gradient, which is generated by the sodium-potassium pump on the surface. In Eq. 7, the strength of 
sodium-potassium pump at surface depends on the ion's concentration, which leads 

a Na —3—, a K = - 2^, a cl = 0, (C.l) 

e e 


where 


Ip — Imax 1 


c 


Na 


c 


Na 


K 


Na 1 


c* 

v --n 


C£ + K k1 


Imax 2 


c 


Na 


t-Na I is 

C,-„ + K Na2 


c* 

v —n 


Q K - 


‘ Kk2 


(C.2) 


We assume that the velocity at surface determines the characteristic velocity scale for the problem. We have 
ion fluxes in the intracellular, extracellular region in Eq. 5 and trans-membrane source of ion in Eq. 6 for 
ion Na + ,K + , Cl~. 

At boundary of the intracellular space, due to the ion pump in Eq. C.l and assumption of conductance at 
surface that G Na — G cl — 0 [51,47], we have 


jNa _ Na jK 
Jin a ' 'in 


js + a K , 


J a — 0 

Jin 


(C.3) 


Since g K — 0 inside of the lens, we obtain 

jf + a K = 0. (C.4) 

This assumption obviously will have to be replaced in applications to other tissues, with a less particular 
distribution of channel proteins. 

By the conservation of fluxes for each ion in Eq. 4, we get 


Ji„ — ^oJex' i = Na,K,Cl, 


(C.5) 
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where <5 q = Therefore, Eq. C.3 becomes 


-5 0 J? x a = a Na , -8 0 jf x = 0, -6 0 J%= 0. (C.6) 

Adding up all three fluxes in Eq. C.6 and since in the extracellular region each ion diffusion coefficient are 
at the same level of approximation, i.e. 

Dg X — O (D ex ), i — Na, K, Cl, (C.7) 


and based on Eq. 10 , we get 


i r n ^ r\ i ^ex^c^Q d ]Wj 

Oex^in T 8oD ex T c -j-Oex H ) V Pe X ~7~ipe X — 0 
dr k B T dr 


(C.8) 


The strength of the ion pump a Na depends on the ion concentration in Eq. C.2 . We choose the scale of a 


Na 


is n Na * based on an experimental estimation [51]. Using Eq. C.8, we take the scale for O m ,, x and Uj n to be 


O* and u* n as 


O* = 2 ( C, Nfl + C 0 K ) , 


„Na* 


11 in q* 


(C.9) 


By mass conservation expressed in Eq. 1, we naturally get the scale of u ex as 


* f—1 * 

u ex ~ ^0 u in- 


(C.10) 


Furthermore, <p * = is used for the scale of electric potential (pj n and cf) ex . For the extracellular velocity 

in Eq. 2, we have 

k _ d ~ k„T d ~ 

(C.U) 


k B T d 


u ex u ex — „,T> r cPex jyPex k e T c l Jf'k ex ' 


}iR 

We think the ^ P ex term balance the velocity u ex . The scale for extracellular pressure P xx is then choose 


P* — 


P Ru *ex 

KexTc 


Therefore, we get 


d ~ d ~ 

U ex — ~F^Pex 8\ -p^fpex/ 

dr /A "' 


dr 


(C.12) 
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where 8 \ — T . For the intracellular velocity, we have 


K in P j n d ~ Ki n fl m k B TO* d 

/<R tfr + /.(R (lr ' r 


(C.13) 


We claim term dr p m and balance at the same level. Therefore, we choose the same scale for the 

intracellular and extracellular pressure, namely. 


p* — p* — p* 

1 1 in L ex' 


ThenEq. C.13 becomes 


d ~ d ~ 

^ 2 11 in ~ ~^ 3 flp P ' n dfOinr 


(C.14) 


where 


<52 — 


V Ru t 


Kin7mk B TO 


:/ ^3 — 


7mk B TO* 


In all, the fluid system Eq. 1 becomes 


Uex — U[ n , 

1 d 


^4y2 fly — $3 (Pex P hi^ + Oex^ r 


(C.15) 


with boundary condition 


where 


Pex — 0 , 

^ 5 U{ n ^3 P in (^in &ex I / 




M;„ut 


:> = 


RM v L m j m k B TO* L s fl s k B TO* 

Appendix C.2. Ions circulation 

The velocity scales and diffusion coefficients in the extracellular and intracellular space are at different 
levels of approximation in our approach. In the following, we put the characteristic diffusion coefficients 
at intracellular and extracellular region and scale of concentration as 


D* t~\CI —- pP pjC/ p* pNa i r-'K 

u ex — u ex T cr — u [ n , ■ 
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In this way, we get Peclet number in the extracellular and intracellular and dimensionless Nernst potential 


as 


p e . — U in R p e _ U ex R gz — — log | — 

£w “ Df ' Pe « - D* ' E - l0§ CL . 


Because g Na — 0 inside of lens, we have K + system as in Mathias's model [51], 


— — (r 2 (p e C K u - D k (—C K ■ 

r 2 dr \ \ £ex ^ exl ex ex l dr^ ex 


-z K C K —d 

ex dr 


(V 2 f p e - r K n- — d k ( —c K + - K c K — &>■ 

r 2 dr \ V etn ^m m l m dr xlt 


= 0, 
-0, 


(C.16) 


with boundary condition 


r K — r K 

ex '-'0 / 

Pe- C K u- -D k ( —C K + z K C K —fi- 

1 t 'in'~in u in ^ '-in j^Y" 


= ^(fc-Eh+2*, 


and Cl system as 


r [f (PeexC%u ex - D% (^C%+z a C%^cp ex 


1 d 
r 2 dr 


x2 ^ Pe- C cl u- — D C1 ( —C cl 

1 c in'~ ln u in | Agin 


7 C/pC/ ~ , h 
z dJ‘ ln 


d 


1 d 


^r 2 dk 1 r 


2 I p e r cl u — r> cl I —C cl 

1 eex'-ex u ex u ex l Ag^ex 


M, 


7 CI 


z Cex d^ e 


( pex £ 


Cl 


(C.17) 


with boundary condition 


c£x — 0^° + Of + ^7 (<pin — <pex) / 

Pe- c a u- - D° (—C. cl +z a C a —d> I -0 

1 e inC- ln U m L- , m I ^~C- m TZ '-'n! ( gg ( r m I U ' 


where 


_ G K k B TR K _ a K R ~ _ M v g a k B TR 2 
s e 2 D* n C* ’ D * n C*' v M ex e 2 D* x C* ’ 


■MexD^x 
M- n* ' 

jv hn u in 


The concentration of Na + can be solved from the following equations 


L^ln + ZTT = 4 (‘Pin ~ <Pex) , 

v i 


Lz'C'ex = ~ S 7 (^in ~ <Pex) / 


(C.18) 
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where 


(C.19) 


n _ M v C m kgT _ MvCrnksT 

6 “ e 2 C*i] ' 7 “ e 2 C*(l — t]) ’ 

From Eq. 11 and use the fact z N,! — z K — 1 and assumption that g Nn = g CI and G Na — G cl 


have 


1 d 
7 2 dr 


Y P^expexUex Dex% j~C ex CFex ^tyex j —(/Pin (pex) 


- E Na - E C1 


< Jl ^ y e inPin u in ^in °~ in J 

„ 1 d f i ~ j : d ^ d 

~ I r ( P e exPexU ex ~ \_,D ex Z j~C e x ~ CT ex ^(p, 


ex / 


with boundary condition 


<pex — 0, 

P e inPin u in — E^E Z in ~ ^in^—tyin — (pin ~ £ ^ + Ip, 


where 


and 


Pin — |z| + S 6 ((pin (pex) / Pex — $7 ((pex (pin ) / Ip — i-.* p* / 

'in (’ 2 'in 


= E D m ( z ’) 2 C in' ^ = E D ex( z ‘) 2 C! 


0, we 


(C.20) 


(C.21) 
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